Non-singular cosmology in a model of non-relativist ic gravity 
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We present a model of non-relativistic gravitational theory which is power-counting renormalizable 
in 3+1 dimensional spacetime. When applied to cosmology, the relativity- violation terms lead to 
a dark radiation component, which can give rise to a bounce if dark radiation possesses negative 
energy density. Additionally, we investigate a cyclic extension of the non-singular cosmology in which 
the universe undergoes contractions and expansions periodically. In both scenarios the background 
theory is well defined at the quantum level. 
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I. INTRODUCTION 



Relativity is commonly assumed to be the foundation 
in constructing models of particle physics and gravity. 
Yet, theories based on relativity often suffer from some 
theoretical defects. Namely, a quantum field theory of 
general relativity cannot be well established since it is 
unable to be renormalized. Various attempts to break- 
ing relativity have been intensively discussed in the lit- 
erature. Recently, motivated by Hofava pj,[2j, models of 
non-relativistic quantum field theory were studied, not 
only theoretically @, H |j| .but also in experimental de- 
tections HSH (see also @|). 

In the present work we are interested in constructing 
a model of non-relativistic gravity and study its cosmo- 
logical implications. This model is power-counting renor- 
malizable in 3+1 dimensions and hence ultraviolet (UV) 
complete. Moreover, its action can recover the exact 
Einstein-Hilbert form in the infrared (IR) limit, and so 
general relativity and Lorentz symmetry in local frame 
do emerge at low energy scales. A generic feature of 
this model is an existence of dark radiation for which 
the energy density can be either positive or negative. In 
the frame of a Friedmann-Robertson- Walker (FRW) uni- 
verse, dark radiation with negative energy density can 
give rise to a bouncing solution [Io| . since it breaks en- 
ergy conditions if accompanied by normal matter com- 
ponents This is the so-called "quintom scenario" 
[H, 0j Il4l Il5l| in which the equation-of-state of the uni- 
verse crosses the cosmological constant boundary. A re- 
markable point of this model is that the violation of en- 
ergy condition does not bring the quantum instabilities 
[l6l[l~7j which often exist in usual quintom scenario. The 
scenario of bouncing cosmology has been investigated 
in models motivated by various approaches to quantum 
gravity [H [H, US HE H3, HI, and analyzed using ef- 
fective field techniques by introducing energy-condition- 
violating matters [ll, 24, 25,(26[, while the generation of 



'Electronic address: 
t Electronic address: 



caiyf@ihep.ac.cn 



msaridak@phys.uoa. gr 



perturbation was studied in [23, [28|, [29|, |30|, |3J, [32j, [33j, [3J 
(and we refer to [351 ] for a comprehensive review). One 
model of bouncing cosmology with a matter-dominated 
contraction was found to be able to provide a scale- 
invariant spectrum |3ol l37l . l38l l39l l4fj| and sizable non- 
Gaussianities jU, |42|. which may be responsible for the 
current cosmological observations. 



To extend, we find that the evolution of a universe in 
a model of non-relativistic gravity might be free of any 
singularities. One example of this scenario is an oscillat- 
ing universe in which the universe experiences a sequence 
of contractions and expansions [43|, 0, |45[ . As shown in 
[44j . for a pivot process (bounce or turnaround) to oc- 
cur, one must require that the Hubble parameter vanishes 
while its time derivative is non- vanishing at that point. 
Thus, it needs to be satisfied that the equation-of-state w 
of the universe is much less than —1 at the bounce point 
while much larger than 1 at turnaround point. Conse- 
quently, both the bounce and turnaround break certain 
energy conditions. An oscillating universe realized in a 
spatially flat universe has bee studied in [44[ by mak- 
ing use of a quintom matter with a quantum instability. 
However, in the current work, we obtain a cyclic scenario 
within a non-relativistic gravity model, which is well de- 
fined at the quantum level but it merely requires that the 
universe is not spatially flat. 



This paper is organized as follows. In Section [H] we 
begin with the quantum theoretical puzzles of usual Ein- 
stein's gravity, and we construct a non-relativistic grav- 
itational model which is power-counting renormalizable. 
In Section IIIII we apply this model in the frame of an 
FRW universe and then we obtain a bounce solution, in 
which we take a canonical scalar field to mimic the nor- 
mal matter component. In addition, we find that the 
cosmological evolution in this framework can be oscil- 
latory and singularity-free, and the corresponding form 
of the scalar-field potential can be straightforwardly re- 
constructed. In Section IIVI we present a discussion on 
the cosmological fluctuations in the model. Finally, in 
section [V] we summarize our results. 
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II. A MODEL OF NON-RELATIVISTIC 
QUANTUM GRAVITY 

We start with a discussion on the main obstacle against 
the usual approach to quantum gravity. In the context 
of quantum field theory all successes are established on 
a solid construction of a perturbatively renormalizable 
model, namely the SU(2) Standard Model, in 3+1 di- 
mensions. However, this procedure does not work in the 
gravity sector, since the gravitational coupling constant 
G possesses a negative dimension ([G] = —2 in mass 
units), which introduces a UV incompletion to the the- 
ory. 

The Einstein-Hilbert gravitational action of is given by 

s g = ^ J dtd 3 x^N (k 13 k^ -K 2 + R^j, (1) 



where 



(2) 



is the extrinsic curvature and R is the three-dimensional 
Ricci scalar. The dynamical variables are the lapse and 
shift functions, N and iVj respectively, and the spatial 
metric gij (roman letters indicate spatial indices), in 
terms of the ADM metric 

ds 2 = -N 2 dt 2 + g tj (dx l + N'dt) (dx> + N j dt), (3) 

where indices are raised and lowered by . 

Attempts on solving the UV incompletion have been 
intensively studied in the literature. Motivated by a Lee- 
Wick model [46|, [47} which shows an improved UV behav- 
ior [48| , gravity involving higher-derivative terms may be 
applied to provide a UV completion. However, this ap- 
proach suffers from the existence of unbounded from be- 
low energy state, and thus its quantization becomes un- 
reliable. Another path to quantum gravity is to build a 
non-local theory [49J , using infinite high-derivative terms 
[HoL l5ll . [521 . [53| , or string theory [54], [55[ , attempts which 
are still in proceeding. 

Motivated by a recent work one realizes that a 
model of power-counting renormalizable gravity may be 
achieved just by adding higher-order spatial derivative 
terms. One important peculiarity in this type of models 
is that Lorentz symmetry has to be given up but it may 
appear as an emergent one at low energy scales. The orig- 
inal model, which is the so-called Hofava gravity, since 
is required to satisfy a 'detailed-balance' condition refer- 
enced from condense matter physics, it still suffers from 
problems such as the over-constraining in UV region, be- 
ing not compatible with current observations even in the 
IR limit. 

Concerning to above issues, the logic of effective field 
theory suggests that a complete action of gravity could 
include all possible terms consistent with the imposed 
symmetries, and the dimensions of these terms ought to 
be bounded due to renormalization. In the frame of 3+1 



dimensional spacetime, a renormalizable term may allow 
6th-order spatial derivatives at most, as pointed out by 
[l[. As a consequence, one could add a modified action 
which involves all the permitted terms: 



A5 



X — J dtd 3 x^/g~N ^> 1 /,•, ,//" + a 2 R 2 



16 

•n,V,/>',; V'/,''' 
+a 5 Vii?V l i? 



n,V,//,;V ! /r 



(4) 



Adding these terms into action (fTJ) and requiring the signs 
of the 6th-order spatial derivatives to be negative, we 
obtain a non-relativistic gravity theory which is power- 
counting renormalizable and the dispersion relation is 
bounded from below. Note that in the following we are 
interested in the cosmological applications of this gravi- 
tational background, that is in its IR limit, and thus we 
do not present the full "running" formalism but we re- 
strict ourselves to the IR expressions. Furthermore, with- 
out loss of generality, we will focus our investigation on 
the quadratic terms, which preserves parity and Poincare 
symmetries, although ones can straightforwardly exam- 
ine a model with all the coefficients included. 

We can insert a matter component in the construction, 
attributed to a canonical scalar field <f> described by 



dtd 3 x^N[-d^d^ - V{cj>)] . (5) 



In addition, we focus on the cosmological frame with an 
FRW metric, 



N= 1 , s«=a 2 (*)7« , JV* = 



with 



dx l dx J = 



dr 2 



! dn 



2 > 



(6) 



(7) 



1 — kr 2 

where k = —1,0,1 correspond to open, flat, and closed 
universe respectively (with dimension [k] = 2 in mass 
units). Both the scalar field and the metric are assumed 
to be homogenous, with their backgrounds being func- 
tions of cosmic time t. By varying N and g^, we obtain 
the Friedmann equations as follows 

8nG r I 

„ Prn + Pk + Pdr 



H 2 = 



H 



-H 2 = -AttG 



1 



1 



Pr, 



;Pk 



'.Pdr 



(8) 
(9) 



where H = — is the Hubble parameter, and the matter 
pressure and energy densities are expressed as 



Pm = \-V{4>) 



Pm 
Pk 
Pdr 



Y + v{4>) 

3k 



8irGa 2 ' 
3(ai + 3a 2 )k 2 
AnGa A 



(10) 

(11) 
(12) 

(13) 
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respectively. In addition, the energy density for mat- 
ter component satisfies the continuity equation p m + 
3H(p m + p m ) = following from the action ([5]), which 
leads to the equation of motion for the scalar field 



ZH<t> + V.4, = , 



(14) 



where the subscript denotes the derivative with re- 
spect to 4>. 

In summary in Friedmann equation |jHJ) we have ob- 
tained the energy density for usual matter component 
(with equation-of-state parameter — 1 < w m < 1, where 
w m = p m / ' Pm) and the usual spatial curvature term (of 
which the equation-of-state parameter is = —1/3 as 
can be seen from ([9])). Besides, we have derived a "dark 
radiation term" , which evolves proportionally to a -4 and 
thus its equation-of-state parameter is Wdr = 1/3 (as can 
be seen from (j^J). This last term reflects a novel feature 
of a model of non-relativistic gravity, leading to smooth 
cosmological evolutions with the initial singularity re- 
placed by a big bounce, if w rn < Wdr and ct\ 4- 3ct2 > 
and the universe is not exactly flat. 



III. COSMOLOGICAL EVOLUTIONS 

In this section we will study in detail the background 
evolution in the context of non-relativistic gravity. In 
particular, we study a non-flat universe with k = ±1 and 

3(ai+3a 2 )fc 2 



4ttG 



we define a critical energy density as p c 
which needs to be larger than zero. Furthermore, we 
normalize a = 1 at the bounce point. 



A. A bounce solution 

Let us take a first look at how it is possible to obtain a 
bouncing cosmology in this model. We assume that the 
universe starts in a contracting phase and that the con- 
tribution of the normal matter dominates over that of the 
dark radiation. This will typically be the case at low en- 
ergy densities and curvatures. As the universe contracts 
and the energy density increases, the relative importance 
of dark radiation compared to the normal matter will 
grow. From ((SJ) it follows that there will be a time when 
H = 0. This is a necessary condition for the realization 
of the bounce. By making use of the continuity equations 
it follows that at the bounce point H > 0. Hence, we can 
acquire a transition from a contracting to an expanding 
phase, which is a cosmological bounce. 

Now, as a specific example, for the scalar field respon- 
sible for matter we consider a mass potential V(<f>) — 
^m 2 (f> 2 . We begin the evolution during the contracting 
phase with a sufficiently large scale factor, so that we ex- 
pect the contribution of the dark radiation to the total 
energy density to be small. Besides we require the con- 
tribution of the spatial curvature term to be negligible 
too (we leave the discussion on this component for the 



next section). For these initial conditions the scalar field 
will be oscillating around its vacuum, and the equation of 
state will hence be that of a matter dominated universe 
with its Hubble parameter being 



V 1 it 



(15) 



on average. In fact, as follows from the Klein-Gordon 
equation l|14p. which in this case reduces to 



tf> + 3H<f> + m 2 4> = , 
one well-known approximated solution is 



0(t) 



1 sin rat 



V3^G mt 



(16) 



(17) 



When t is smaller than m -1 , the scalar field cf> is oscil- 
lating and the universe behaves as a matter-dominated 
one. Once t ~ m _1 , the amplitude of <f> reaches the 
Planck scale and the scalar enters a "slow-climb" phase. 
During this period, the energy density of the scalar varies 
very slowly but the scale factor decreases rapidly. There- 
fore, the dark radiation will become similar to the scalar 
field density very soon. Once there is p m — p cr , we 
will obtain a big bounce. Since in this case the uni- 
verse has experienced a matter-dominated contraction, 
this scenario is the so-called "matter-bounce" . 



B. Cyclic scenario 

Having investigated the realization of bouncing cos- 
mology in the present non-relativistic gravitational 
framework, we now study a sub-class of cosmological evo- 
lution without any singularities, that is a realistic and 
physical cyclic scenario. 

A cyclic scenario could be straightforwardly obtained 
if we consider a negative dark radiation term and a neg- 
ative curvature one. According to t\l'2\i and (|13p . this 
is achieved if ot\ + 3«2 > and if k = 1, that is for a 
closed universe. For simplicity we consider the matter 
component of the universe to be dust, that is to pos- 
sess w m — 0. During the expansion, the energy densities 
of all components are decreasing. However, the energy 
density of the curvature term decreases much slower rel- 
atively to others. Thus, its contribution will counterbal- 
ance that of dark matter, triggering a turnaround with 
H = and H < 0, after which the universe enters in the 
contracting phase. As described in the previous subsec- 
tion, after a contraction to sufficiently small scale factors 
the dark radiation term will lead the universe to experi- 
ence a bounce. Therefore, the universe in such a model 
indeed presents a cyclic behavior, with a bounce and a 
turnaround at each cycle. 

A more careful analysis reveals that the negative cur- 
vature term, that is a closed universe with k = 1 , is not a 
necessary condition. Indeed, its role on the competition 
of the positive dark matter density can be equivalently 
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fulfilled by a small negative cosmological constant added 
to the potential V(<f>). In other words, a slightly negative 
V{(j)) can trigger the turnaround at large scale factors, 
even if the universe is open, i.e with k = — 1. 

The above general examples reveal the realization of 
a singularity-free cyclic universe in a framework of non- 
relativistic gravity. This behavior can be obtained for 
both closed and open universe, and the only restriction 
is the requirement of a matter content with equation- 
of-state parameter between — i and i. Thus, having 
extracted its general features, in the following we will 
explicitly construct the class of models that allow for 
cyclicity, and in particular we wish to appropriately re- 
construct the corresponding scalar potential. 

Let us first start from the desired result, that is to im- 
pose a known scale factor a(t) possessing an oscillatory 
behavior. In this case both H(t) and H(t) are straight- 
forwardly known. Therefore, we can use the the Fried- 
mann equations ©,(© together with (flU )) . (fTTj) . in order 
to extract the relations for <fi(t) (through 4>(t)) and V(t), 
acquiring: 



0(4) 
V(t) 



± 



dt'\ 



H(t') 



H(t) 
8irG 



4ttG 3 
3H 2 {t) 2 



Pk{t') - -PdrW), (18) 



8vrG 



Pk(t) - -Pdr(t) 



(19) 



Note that the a(4)-form or the parameter-choices must 
lead to a positive 4> 2 {t). Finally, eliminating time be- 
tween these two expressions we extract the explicit form 
of the potential V((f)). Thus, performing the procedure 
inversely we conclude that this particular V(4>) generates 
the desired a(4)-form. 

We now proceed to a specific, simple, but quite general 
example. We assume a cyclic universe with an oscillatory 
scale factor of the form 



a(t) = Asm(cut) + a Cl 



(20) 



where we have shifted t in order to eliminate a possible 
additional parameter standing for the phase. Further- 
more, the non-zero constant a c is inserted in order to 
eliminate any possible singularities from the model. In 
such a scenario t varies between — oo and +00, and t = 
is just a specific moment without any particular phys- 
ical meaning. Finally, note that the bounce occurs at 
^b(4) — o- c — A, which can be set to 1. Straightforwardly 
we find: 



H(t) 



H{t) 



Aui cos(wi) 
A sm{ujt) + a c 
Ato 2 [A + a c sin(w4)] 
L4sin(w4) + a c ] 2 



(21) 
(22) 
gives 



and thus substitution into (12]) , (13]) , UB and vM 
the corresponding expressions for (jj(t) and V(t). 

In order to provide a more transparent picture of the 
obtained cosmological behavior, in Fig. Q] we present the 



evolution of the scale factor (|20]) and of the Hubble pa- 
rameter (|2"Tj) with A = 10, u> = 0.1 and a c = 11, where 
all quantities are measured in units with 8irG — 1. This 
choice is consistent with our setting a = 1 at the bounce. 
In Fig. [2] we depict the corresponding behavior of (j>(t) 
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FIG. 1: The evolution of the scale factor a(t) and of the Hub- 
ble parameter H(t) of the ansatz H2(J[I , with A — 10, u> = 0.1 
and a c = 11. All quantities are measured in units where 
8irG = 1. 



and V{t) for the scale factor of Fig. [T] in the case of a 
closed universe (k = 1), and for model parameters ax = 1 
and cti = 1 (in units where 8ttG =1). Finally, climinat- 
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FIG. 2: (fr(t) and V(t) for the cosmological evolution of Fig.\7\ 
in the case of a closed universe (k — 1) with cti = 1 and 
ot2 = l. All quantities are measured in units where SixG = 1. 
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ing time between <fi(t) and V(t) allows us to re-construct 
the corresponding relation for V(<f>), shown in Fig. [31 



V 




FIG. 3: V((j>) for the cosmological evolution of Figs. Q] andfJI 
All quantities are measured in units where 8tyG = 1. 



In this case <f> = sin 1 ( V ^y^ Vc J , where V((fi(t)) = 

V(t) = Q2(a,a,a) with Q2(a,d,a) the right hand side of 
expression (TiUl) . Therefore, differentiation leads to: 



i 



V u>v 



Q-2{a,a,a) — V c 
Vo 



2 dt 



i(a,a,a)] (25) 



and thus we obtain 



Qi(a,d,'d) 



V U! V 



2 (a,a,a) — V c 



Vo 



2 dt 



[Q 2 (a,a,a) 



where as we have mentioned, Qi(a, d, a) is the right hand 
side of expression (fl8|) . 

Differential equation ([2"o]) cannot be handled analyti- 
cally, but it can be easily solved numerically. In Fig. [4] we 
depict the corresponding solution for a(t) (and thus for 
H(t)) under the ansatz (JUj) with V = 5.25, uj v = 0.25 
and V c = 5.25, with k = 1, ax = 1 and 012 = 1 (in units 
where 8irG = 1). The potential parameters have been 
chosen in order to acquire a cyclic universe with a(t) ~ 1 
at the bounce. 



From these figures we observe that an oscillating and 
singularity- free scale factor, can be generated by an oscil- 
latory form of the scalar potential V(4>) (although of not 
a simple function as that of a(t), as can be seen by the 
slightly different form of V(4>) in its minima and its max- 
ima). This V A (</))-form was more or less theoretically ex- 
pected, since a non-oscillatory V(<p) would be physically 
impossible to generate an infinitely oscillating scale fac- 
tor and a universe with a form of time-symmetry. Finally, 
we stress that although we have presented the above spe- 
cific simple example, we can straightforwardly perform 
the described procedure imposing an arbitrary oscillat- 
ing ansatz for the scale factor. 

The aforementioned bottom to top approach was en- 
lightening about the form of the scalar potential that 
leads to a cyclic cosmological behavior. Therefore, one 
can perform the above procedure the other way around, 
starting from a specific oscillatory V((f>) and resulting to 
an oscillatory a(t). In particular, (jTSJ) is written in a 
compact form as </» 2 (i) = Qi(a, d, a) and similarly (|19p as 
V(t) = Q2(a,d 1 a). Thus, we can invert the known form 
of V(4>) = V(0(tj) obtaining <p(t) = (Q 2 (a,d,d)). 

Therefore, <j) 2 (t) = {j- t [V^ (Q 2 (a, d, a))] } 2 . In con- 
clusion, the scale factor arises as a solution of the differ- 
ential equation 



Qi(a,d,d) = \ - [V { - 1} (Q 2 (a,d,d)) 



(23) 



As a specific example we consider the simple case 

V(4>) = V sin(wy 0) + V c . (24) 




FIG. 4: The evolution of the scale factor a(t) and of the Hub- 
ble parameter H(t), for a scalar potential of the ansatz \2J$ 
with Vo ~ 5.25, u>v = 0.25 and V c = 5.25, with k = 1, 
ax = 1 and ct2 = 1. All quantities are measured in units 
where 8nG = 1. 



Let us now present the corresponding simple cyclic ex- 
ample in the case of an open universe (k = — 1). In Fig.[S] 
we depict the evolutions of the scale factor a(t) and thus 
of the Hubble parameter H(t), under the ansatz (f2"4"]l with 
Vo = 3.15, lo v = 0.25 and V c = 3.13, with fc = -1, a x = 1 
and a 2 = 1 (in units where 8irG =1). As described 
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this issue in future studies. 
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FIG. 5: The evolution of the scale factor a(t) and of the Hub- 
ble parameter H(t), for a scalar potential of the ansatz 
wtth Vo = 3.15, uy = 0.25 and V c = 3.13, with k = — 1, 
ax = 1 and ct2 = 1. All quantities are measured in units 
where SttG = 1. 



above, in the case of an open universe one needs the scalar 
potential to be negative for field values corresponding to 
large scale factors, in order for the turnaround to be trig- 
gered. However, since at that regime the curvature term 
is very small, even very small negative potential values 
can fulfill this condition, as can be seen by the specific 
example of Fig. where Vq = 3.15 and V c = 3.13 leading 
the minimal value of the potential to be —0.02. 

Note that the specific examples of the above figures are 
just simple representatives of cyclic behavior in our grav- 
itational and cosmological construction, and they corre- 
spond only to a sub-class of the whole set of cyclic mod- 
els. Obviously, one can straightforwardly generalize the 
aforementioned procedure in any periodic model. For 
instance, imposing any periodic function for a(t) one re- 
sults in the corresponding periodic scalar potential V (</)). 
Similarly, imposing any periodic potential V(<f>) one can 
solve the differential equation (f2"3"| and extract the result- 
ing periodic a(t). 

We close this section by mentioning that, although the 
bounce solutions arise owing to the presence of a dark 
radiation component with negative energy density, they 
can also be obtained if ordinary radiation with positive 
energy density is present. When ordinary radiation is 
involved, it has to be generated from the reheating pro- 
cess of a primordial field namely inflaton or <f> appeared 
in our model. Therefore its domination only takes place 
after reheating of which the energy scale is much lower 
than the bounce scale. Moreover, in the late time evo- 
lution normal radiation would be erased during matter 
dominated period, and hence will not affect the bounce 
solution in next cycle. We will address on the details of 



IV. FLUCTUATIONS THROUGH THE 
BOUNCE 

A model of non-relativistic gravity is usually able to re- 
cover Einstein's general relativity as an emergent theory 
at low energy scales. Therefore, the cosmological fluctu- 
ations generated in this model should be consistent with 
those obtained in standard perturbation theory in the IR 
limit [loj . This result has been intensively discussed in 
the literature (see e.g. [Hj]). In particular, the pertur- 
bation spectrum presents a scale-invariant profile if the 
universe has undergone a matter-dominated contracting 
phase [26|, [4(J E2] • However, the non-relativistic correc- 
tions in the action (H|) could lead to a modification of the 
dispersion relations of perturbations. This issue has been 
addressed in [56j 1 , which shows that the spectrum in the 
UV regime may have a red tilt in a bouncing universe. 
Moreover, the perturbation modes cannot even enter the 
UV regime in the scenario of matter-bounce. So the anal- 
ysis of the cosmological perturbations in the IR regime 
is quite reliable. 

Things become complicated but more interesting in a 
cyclic scenario. Usually, a particular perturbation mode 
in the contracting phase is dominated by its growing ten- 
dency, but in the expanding stage it becomes nearly con- 
stant on super-Hubble scales. Therefore, the metric per- 
turbation is amplified on super-Hubble scales cycle by 
cycl e 15911 , and also the slope of its spectral index is vary- 
ing [6fJ. However, it is known that the contribution of 
fluctuations has to be much less than the background 
energy. This prohibits the metric perturbations to enter 
the next cycle if 8p/ p ~ O(l), unless the universe can be 
separated into many parts independent of one another, 
each of which corresponding to a new universe and evolv- 
ing up to next cycle, then separate again and so on. In 
this case, the model of cyclic universe may be viewed as 
a realization of the multiverse scenario [H, [6l], [62j • 



V. CONCLUSIONS 

In this work, we have studied the possibility of con- 
structing a model of non-relativistic quantum gravity in 
3+1 dimensional spacetime. The novel features of the 
gravitational sector are reflected in new terms that are 
present in the IR, that is in the cosmologically interesting 
regime. Our results show that this model can give rise to 
a non-singular cosmology, for which the initial singular- 
ity is replaced by a big bounce if we require that the dark 



1 we refer to Refs. [571. l58|| and references therein for the perturba- 
tions of a pure expanding universe in Hofava-Lifshitz cosmology. 
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radiation term is negative. Specifically, we have consid- 
ered an example in which the normal matter component 
is a free scalar field. In this case we have obtained a 
bouncing universe with a matter-dominated contracting 
phase, and so it may be responsible for the formation 
of a scale-invariant primordial spectrum. To extend, we 
have also investigated the realization of a cyclic scenario, 
by re-constructing the potential of the scalar field which 
leads to a universe with an oscillatory scale factor. 

Recently, the scenario of oscillating universe (originally 
proposed by [63[ and awaked later as ekpyrotic/cyclic by 
(64U65', 66]), in which the universe experiences a sequence 
of contractions and expansions, has been widely studied 
in the literature, for instance in the context of loop quan- 
tum gravity [i3,E3,[sl,[S§|j including matter components 
violating energy conditions 7(1 7l|, [?2|, [iD, [7J, [75| , in the 
frame of strin g co smology |64U76ll and within the brane- 
world 0, [zllzl (see Refs. [M M H H M, HE 
HH, H3 for recent developments on oscillating universes). 
The main peculiarity of the current work is that the os- 
cillation depends on the dark radiation term, which is 
present only for not spatially flat geometry. However, 
the main advantage is that the background theory is well 
defined at the quantum level. 

In conclusion, we see that the present model of power- 
counting renormalizable, non-relativistic gravitational 
theory in 3+1 dimensional spacetime, can naturally lead 
to a bounce and to cyclicity as particular sub-classes of its 
possible induced cosmological behaviors. The fact that 
the theory is UV complete and that instabilities do not 
arise at the quantum level, makes future investigations 
on its cosmological implications quite interesting. 

As an end, we would like to comment on a principal 
difference between our model and the intensely studied 
Hofava-Lifshitz cosmology [H H|. In a model of Hofava- 
Lifshitz cosmology (with or without detailed balance con- 



ditions), it is claimed by [88| that, there is one extra 
scalar modes in its perturbation theory which becomes 
strongly coupled as the parameters approach a desired IR 
fixed point. To understand this point, we would like to 
recall that such a strong coupling problem only exists in a 
system of which the original scalar modes have obtained 
effective mass terms and so lead to the generation of an 
extra degree of freedom, as what we have understood 
in quantum field theory very well. As a consequence, 
modifications on general relativity often suffers from this 
problem, namely, in the theory of Pauli-Fierz massive 
gravity [89| the longitudinal scalar becomes strongly cou- 
pled when the mass approaches zero [90| . leading to the 
famous vDVZ discontinuity [9l|, Hofava-Lifshitz 
cosmology also suffers from the strong coupling prob- 
lem, since the theory manifestly contains parity- violating 
termsfHI, [94[ which bring effective masses for gravitons. 
However, this does not happen in our model since the 
spatial curvature terms and its spatial derivatives intro- 
duced phcnomenologically only contribute to the disper- 
sion relations of two polarizations of gravitons without 
bringing effective mass terms. In this case, there is no ex- 
tra degree of freedom in our model and the longitudinal 
part of metric perturbations can be fixed by Hamiltonian 
constraint when combined with matter component. 
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